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SUMMARY 

given to the resulting el fectiveness ol aeolian erosion. 

INTRODUCTION 

conditions tor initiating particle movement by the wind with discussions based mainly on labora 
tory simulations. 

Wo wilt 1 St brietly review the current understanding of martian dust storms, then describe the 
We will I..st uiuiy review I K Research Center to simulate martian aeolian 
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MARTIAN DUST STORMS 


More llum one luiiulreil years of observations o1’ Mars Itave produced an impressive catalop 
documenting the changing patterns on that planet. In recent years stall’ at several observatories have 
closely followed the initiation and growth of global martian dust storms. Summari/.ing from 
reference I , these observations show that major dust storms typically originate in three general areas 
of Mars Hellespont us, Noachis, and Solis Planum all elevated plateaus that are between 20“ and 
40° S latitude (fig. I ); Viking results are showing that local dust storms can occur in many other 
areas as well. Major storms typically begin near or slightly before the time ot southern hemisphere 
solstice (close to perihelion), at the start of martian southern hemisphere summer, or late spring; 
however, some years there are no major dust storms. 

The major storms appear to go through three phases. In phase I, lasting about five days, the 
storms begin as bright spots or cores, about 400 km or smaller in diameter. Phase II is the expansion 
of the storm which can last from about 35 to 70 days. Expansion takes place by having secondary 
cores develop around the primary cores of the first phase, until eventually the entire planet is 
affected. Moreover, the development of the storm path does not seem to be topographically related; 
for the largest storms, the entire planet may be totally obscured. Phase III marks the decay of the 
storm and lasts from 50 to 100 days. The first areas to clear are the poles and topographically high 
regions, such as the summits of the shield volcanoes. Although major dust storms do not occur 
every year (phase II - expansion - may not fully develop), the occurrence is fairly frequent 
(ref. 1). 

The Mariner 9 mission to Mars in 1971 provided the unique opportunity to observe dust storm 
activity at close range Even before the spacecraft went into orbit. Earth-based observations showed 
that a major dust storm totally obscured the planet. Photographs obtained from Mariner 9 coupled 
with other data produced a wealth of new information on martian dust storms. Average particle size 
in the atmosphere was found to be less than 2 fjtm, or about the same as the particles carried over 
the Atlantic by major Saharan dust storms. Tlie dust was found to be well mixed in the atmosphere 
to heights of 30 to 40 km and had the effect of raising the atmospheric temperature by as much as 
50 K. 

After the 1971 dust storm had essentially cleared, Mariner 9 photographed numerous local 
dust storms and a host of features that are attributed to aeolian processes (ref. 2 and others). Major 
findings include dune fields, polar laminated terrain (possibly composed of alternating iayers of 
windblown dust and dust-ice mixtures) and yardangs. The most abundant type of aeolian feature, 
however, is the crater-streak (fig. 2), of which there are two types - dark streaks and light streaks. 
Dark streaks were seen to appear and change size and shape in a matter of weeks, light streaks 
remain relatively stable, taking years to change (ref. 3). Most dark streaks appear to represent areas 
where aeolian erosion has removed material, while some light streaks appear to represent dust 
deposits. The morphology of many streaks can be explained by the geometry of the wind flow field 
around and over craters and the resulting zones of wind-scour and deposition (refs. 4 8). 

The Viking mission currently in progress is yielding even more data on martian aeolian 
processes. Two Viking Orbiters are producing photographs of the surface over large areas at 
resolutions substantially improved over Mariner 9, while two landers are obtaining the first pictures 
of the surface, making measurements of the wind speeds and directions, and obtaining data on the 















ORIGINAL PAGB IS 
OP POOR QtJAIJTY 

atmospheric dust content for the two different landing sites in the northern hemisphere. Under 
pictures at both sites show accumulation of windblown sediment as dune-forms (fig. as well as 
evidence ol aeolian erosion in the form of possible ventifacts (refs, and 10). Orbiier |)ictures show 
newly discovered dune fields (fig, 4) and numerous crater streaks (ref. II). Aeolian features 






7 — 

7 • ■' - 






m *"'4' y 




(a) Panoramic view of the surface of Mars from Viking Lander I in Chryse Planitia (see fig. 1 ). This image was taken 
in eariy morning (about 7:30 iocal time), August 3, 1977, and covers the view from the northeast at the teft, to the 
southeast at the right. The large boulder at the left is about I by 3 m and is about 8 m from the spacecraft. The 
Ught-colored materials are deposits of windblown particles. The shape of the deposits suggest erosion by winds 
blowing from the upper left to the lower right (ref. 9), 
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(b) View of martian surface from the Viking Lander 2 site in Utopia Planitia (see fig. 1), showing the rock -littered 
terrain and accumulations of fine-grained, windblown material. The horizon, which is some 3 km distant, is actually 
level; the tilt in the picture retlects the 8“ tilt of the spacecraft and the manner in which the panoramic pictures are 
obtained (ref. 10). 


Figure 3. V.ews of martian surface from Viking. 
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Figure 4.- Viking Orbiter 1 view of a sand dune field contained within Gatigis Chasma, a canyon of Vallis Marineris 
(see fig. 1). The field is about 50 km long and up to 20 km wide. Individual dunes are about 500 m across. The 
illumination is from the right toward the left. 

observed from both orbit and on the ground will be monitored throughout the extended Viking 
mission into 1978 for observation and for correlation with the meteorology experiment results 
(refs. 12 and 13) and other Viking data. For example, local dust storms have been observed 
developing in the southern hemisphere (figs. 5 and 6) earlier than expected. By late martian spring 
(1977), most of the southern hemisphere was obscured, although the southern polar cap remained 
clear. 

From Earth-based observations and numerous spacecraft missions, it is apparent that aeolian 
processes have played an important role in altering the surface of Mars. Knowledge of the physical 
and to some extent the chemical conditions that govern aeolian processes on Mars is paramount to 
our understanding of both the present environment of the planet and of its geological history. 

In order to understand the nature of aeolian processes on Mars, a collaborative research 
program was established in 1972 that is a multidisciplinary effort involving a geologist, a planetary 
physicist, and aerodynamieists; it includes laboratory simulations (primarily wind-tunnel work), 
field studies of aeolian processes on Earth, and analyses of spacecraft data. Previous results Irom 
this program are reported in references 4 b and 14 20. 
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Figure 5 - Viking Orbiter 2 photomosaic of dust storm (arrow) more than 300 km across 
inside the Aigyre Basin (see Hg. 1 , Argyre Planitla). 
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Figure 6. Viking Orbiter views “before” (left) and “after” (right) of the region west of Vallis Marine.is. The mosaic 
at left was taken July 31 , 1976, and shows the surface in sharp detail; the picture at riglrt was taken March 25, 
1977, and shows the same area covered by diffuse clouds of water ice and dust. 



TUli MARTIAN SURFACIi WIND TUNNI-L (MARSWH ) 


When the possibility of dust storms on Mars was first reali/.ed, many investiwatois attempted to 
predict tile wind speeds and other conditions needed for martian winds to move fine-nrained 
material, In most cases these investigators turned to the classic work of Bagnold (ref. 21) for the 
expressions describing aeolian processes and used the appropriate parameters for Mars in fornutlat- 
ing tlieir predictions. Assuming rouglily analogous geological conditions t)f particle mineroiogy . etc., 
the greatest problem in extrapolations from Ivarth to Mars are differences in temperalure and 
atmospheric surface pressure and, to a tnuch lesser extent, the two-thirds less gravity on Mars, The 
atmospheric surface pressure on Mars ranges from about 1 mb (100 Pa) on mountain tops to about 
10 mb (1000 Pa) -- less than 1/1 00th that of the liarth -• with u nominal average of about 53 mb 
(530 Pa). Titus, the substantially “thinner” martian air, composed of 95 percent CO^. 1 3 percent 
At, 2-3 percent Nj, 0.1 -0.4 percent Oj (ref. 22), must be blowing considerably faster than on 
Earth to accomplish the same effect, or to initiate particle movement. Unfortunately, imprecise 
knowledge of the effect of winds upon a loose particulate surface at low pressure makes extrapola- 
tions to Mars questionable. Widely divergent results for predictions of threshold wind speeds on 
Mars among various investigators prompted another approach, namely, experiments conducted at 
low pressure. 

A low-pressure chamber at NASA-Ames Researcl, Center, constructed in the early 1900’s for a 
series of acoustic and structural tests at low p assure on rockets and since then vacated, was 
recommissioned to house the Martian Surface Wind Tunnel (MARSWtT).. The chamber is a 
pentagon-shaped, concrete tower 30 m high, with a floor area of 164 m^ and a total chamber 
volume of 4058 m^ (fig. 7). Tlie entire chamber can be evacuated to a minimum pressure of T8 mb 
(380 Pa) in a period of about 45 min using a steam-ejection system. A 7.6-m X 7.9-m door permits 
large experimental apparatus to be placed inside the chamber, and numerous plumbing and 
electrical fittings allow a variety of experiments to be conducted at low pressure. MARSWIT 
occupies the center of the chamber, as shown in figure 7, 

MARSWIT (fig. 8) has an overall length of 14 ni, with a 1.1 -m^ test section located 5 m from 
the entrance. The tunnel walls are constructed of 2.4-cni clear Plexiglas to enable ready viewing. 
The tunnel is driven by a network-ejector system consisting of 72 equally spaced 1 .6-m nozzles 
located in the diffusor section. High pressure air (up to 9.86 kg/cm^) is forced througti the nozzles 
to induce flow of air through the tunnel. The maximum attainable free-stream airspeed is 13 m/sec 
at atmospheric pressure, increasing to 180 m/sec at 5 mb (500 Pa). 

A naturally turbulent boundary layer occurs inside the tunnel at atmospheric pressures; 
however, at low pressures (corresponding to the range of martian surface pressures), it is necessary 
to “trip” the boundary layer in the entrance area to ensure that the boundary layer is turbulent. 
Tiiis is accomplished by means of vortex generators, consisting of small pebbles that are secured to 
the tunnel floor. The bed of pebbles extends approximately 3.5 m downstream from the en ranee 
section. The wind tunnel boundary layer that is achieved corresponds to a neutrally stratified 
atmosphere in which the Monin-Obukhov stability length is infinite, hence the ratio of local surface 
roughness height to the stability length is zero. A finite value of the stability length could be 
achieved in the tunnel by heating or cooling the entire tunnel floor to obtain unstable or stable 
stratification, respectively. 
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(a) Schematic diagram of the Martian Surface Wind 
Tunnel (MARSWIT), showing the low pressure 
chamber, the tunnel, and the control room. 


(b) Photograph ol MARSWll' removed from the low 
pressure chamber for maintenance. Right side of the 
tunnel is the diffuser, exit end of the 13-m open 
circuit tunnel. 


Figure 7. - Martian Surface Wind Tunnel (MARSWIT). 
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Figure 8. - Diagram of MARSWIT, showing cntraiiec cone and flow straighteners (Section A-A), Test faction with 
various pressure probes (Section B-B), and diffuser-drive system (Section C-C). 
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Wiiuls Uiinspoil piirtick's by smtaa tiaclion, sus|U'nsiun. ami sallation, tliv latU'i a li'iiii iih’iI 
to dfscribi* homicint? t^iains (dcriviul Imm tlu* l.atin vfil> MiUm\ nraniiit-' “Mo leap or ilaiuv”), 
Mulciial in tiiution imrvcs very close to or on the snrraee by i(tllin(>, sliding, ele., whereas partieles 
in suspension are lu'ld alol't by eddies of turbulent flow. I’ariieles im l-.arth are most easily moved by 
the wind by saltation, with both surface creep and suspension resulting primarily from the inipaet 
of saltatinn ttrains. Thus most research, including im>st of that of Hattnold (ref. - I), has concen- 
trated on the conditions of wind speed, surface roujrhttess, etc., needed to initiate sallation. Hatmohl 
defined two types of threshold static tlireshold and dynamic threshold. Static threshold involves 
wind conditions needed to initiate grain movement for a bed of (rartieles at rest. Dynamic threshold 
involves the initiation of grain movement for a bed that is being impacted by grains saltating from 
upwind; typically, somewhat lower wind speeds are needed for dynamic threshold tlian for sialic 
threshold. The tlireshold experiments described here arc concerncil with static threshold only. 


The basic physics of windblown material is considered to be essentially the same whether it 
occurs on the surface of barth or Mars. The primary difference is due to the comparatively thin 
martian atmosphere which requires much higher wind velocities to accomplish the same effect as on 
Karth. As the wind speed increases over a given surface, there is a minimum wind velocity (r) that is 
needed to initiate motion for particles of a given size and density. Corresponding to this miniimim 
wind velocity is a threshold surface stress (r^) that is needed to set particles into motion. From this 
threshold stress, t^. a threshold friction velocity may be defined as where p is tile 

fluid density. lixperiinentally can be determined from knowledge of the vertical velocity 
structures above the surface. 


Based on an evaluation of the important terrestrial parameters influencing the threshold 
friction velocity, Bagnold (ref. 21) derived the following expression for the tiueshoki friction 
velocity on Farth 


I Kp„ - p)aD„ 

(I) 

where is particle density, g is acceleration due to gravity. is the mean tiarticle diameter, and A 
is an i xperimental eoeffieient. F'or particle motion to occur, the force caused by the horizontal 
wind stress must exceed the gravity force exerted on the particle. In the case of partieles having 
densities much greater than the donsitv of the fluid, equation (1 ) .simplifies to 

A ■■■ ^ 


I'he dimcrisiunless quantity A appearing in these equations was assumed by Bagnold to be a 
unique function of only the particle friction Heynold.s number H. li is defined as 
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wlu’iv V is the kinvnuilU' viscosity, I'lic piiiitiiicicr H proviilcs a mcasmv of how tiirhulcnl the flctw is 
urotMul the particle. In luiiuiple (he detemiinalion of the fiinctioii Adi) is all that is necessary to 
obtain terrestrial values ot if the other piirameters are known (e.a.. Djy p.x, an., r \v<jii!cl 
ordinarily he specified). Ihis function appears to pretlict accurately the values of 1'*^ for I'.arlh as 
lont? as //!> I: however, it may not predict correct resnlts for pressures Uiwer than atmospheric, 
such as those on Mars. 


Thus it is jicnerally possible to predict the values of for a range of conditions on liarth. 
and, for at least a first-order approximation, to |)iedicl for Mars by using the appropriate 
martian values for the parameters involved. Accurate predictions, however, for liarth ami especially 
for Mars involve several complicating factors. Ihe original expressions for particle movement 
derived by Bagnold and amplified by others do not take into account such factors as lift ffcc^, 
interparticle forces due to cohesion, electrostatic charges, etc., or the Magnus effect (a lift functioii) 
from spinning grains. Because these are important not only for understanding aeolian processes j't 
Farth but are essential for extrapolations to Mars, some ol' these factors have been systematically 
investigated in the laboratory for one-atmosphere conditions and reported on earlier frefs. 4 ^ 
and I K 20). 


Despite the refinements that have been made in expressions describing wind threshold 
conditions on liarth, many fundamental c|Uestions remain regarding the behavior of windblown 
material under the low atmospheric pressure conditions on Mars; it was for this reason thai 
MARSWIT was established. Tlie first results from MARSWIT are described in a preliminary report 
by Greeley et al. (ref. 15): in this section, we describe more fully the experimental procedure and 
the threshold results for low pressure conditions. 

There is a basic difference between the laboratory simulations and martian conditions. In the 
wind lunncl the working lluid is “liarth” air; however, on Mars tile atmospheric gas is 95 percent 
carbon dioxide with lesser amounts of argon, nitrogen, and oxygen. The universal gas constant for 
air is /? = 0.287 KJ/KgK, while for CO 2 it is 0.189 KJ/KgK. Hence, for air at a pressure of 5 mb 
(500 Pa) the gas density is b.20 X lO"*’ Kg/m ’. and for (’O 2 at the same pressure the density is 
9,42 X I IT ’ Kg/m’ or 52 percent more than air. Both densities were calculated for a constant 
teniperaturc of 285 K. This represents a substantially larger density of tlie gas on the surface of 
Mars than comparable wind-tunnel tests have at the same pressures and temperatures. This will 
affect the forces on individual grains since the density of gas is directly related to the force. 1'he 
absolute viscosity of air is approximately oiie-fifth greater than that of CO 2 at similar temperatures 
which also complicates the simulation. Thus MARSWIT r.’sults cannot be compared directly to 
similar martian surface conditions, and extrapolations must he developed to take into account these 
differences in order to make predictions for a wide range of surface pressures and temperatures on 
Mars. 
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Tliroshoid friclioi) velocities are typically plotted as a function of particle diameter 
(fig. *)). Several different materials having a range of particle diameters and densities were tested in 
the MARSWIT. Hach experiment consisted of placing a patch of particles in the tunnel test section, 
evacuating the chamber to its minimum pressure (-^3.8 mb, or 380 Pa), and then increasing the 
wind speed through the tunnel until saltation was initiated. The free-strearn wind speed was noted, 
correlated with a standard velocity profile for the tunnel at the given pressure and converted to a 
value of (fig. 10). The chamber pressure was then increased and the procedure repeated for the 
higher value of surface pressure. This process was repeated for each particle test case, and threshold 
velocities were obtained as a function of pressures ranging from ~3.8 mb (380 Pa) to ~l000mb 
ilQO Pa). 
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Figure Comparison of threshold fridion speed vs particle diameter for Farth, Mars, Venus, and in liquid. In all 
cases except liquid, tlicrc is an “upturn'' in Uic speed needed to set smaller grains in motion; thus, on each planet, 
there ir. an opduunn particle size for movement by lowest winds (ref. 36 ). 
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(a) Several velocities are run at each pressure. 
Heiglit of boundary layer is determined. 


(b) Data is nondimensionalizcd and compared to 
standard profiles to insure that it is turbulent. 
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(c) Log slope is obtained to determine K*, from (d) The values of are checked to ensure 

relationship “(f'loo “ V'loVS.TS. linearity. 


Figure 10.- Determination of friction velocity, from boundary-layer proftlcs. 


Test materials - Three kinds of materials were used in the tlireshoid friction velocity experi- 
ments: silica microspheres, natural silt, and finely ground walnut shells. These materials (fig. II) 
were selected because they afforded a variety of particle diameters and densities. 


With knowledge of the geological environment of Mars derived from Mariner data and Viking 
Orbital information plus results obtained from the Viking Lander (refs. 0, 10, 23, and 24). 
windblown particles probably include weathered basalt fragments, clay materials, durientst (crust- 
like material that may consist of agglomerates of fine particles), and perhaps certain “heavy” 
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Figure 11. Hiotomicrographs of some types of the particles used in threshold tests performed under low 
atmospheric pressures. Except for the losses, several size ranges of each type were run in the tests to provide a 
wide range of particle diameters and densities. 

minerals such as magnetite (refs. 25 and 26). Thus the specific gravities of martian windblown 
particles probably range from less than 2 to more than 3. Gravity on Mars is 0.38 that of Barth; 
thus, less force is required to initiate particle lift and subsequent movement. Because gravity is a 
parameter that cannot he controlled in MARSWIT, some other means must be taken to simulate the 
lower gravity of Mars. For threshold simulations conducted on F.arth, particles 0.38 as dense as 
those expected on Mars can be used to offset partly the difference in gravity. The ground walnut 
shells are considered appropriate since their density is about one-third that of possible windblown 
martian material and their shape is similar to natural windblown minerals (fig. 12). Thus the walnut 
shell material was used to derive most of the threshold velocity data tor Mars; silica microspheres, 
silt, and other materials were used to provide a wide range of particle densities for derivation of “A” 
vs “B” relationships. 


Cohesion resulting from water adsorbed on particles raises the threshold velocities above that 
for “dry” particles. To reduce the effect ol‘ this parameter, the test materials were heated in an oven 
to a temperature of about 425 K several hours prior to each run to remove most ot the moisture. In 
addition, the test bed of the tunnel contains a heated plate so moisture would not resorb on the 
particles during pumpdown. The temperature of both the particles and the bed, however, was not 
sufficient to significantly alter the characteristics of the wind boundary layer, hence the Monin- 
Obukhov Length is intlnite and the boundary layer may be eonsidered neutral. 


Another important parameter ol the test particles is the distribution ol particle sizes within 
each sample. Although only limited experiments with mixed particle size samples have been 
conducted, it is noted that the presence of a relatively few la.ge grains can cause threshold at a 
nuicli h)wer value than if the sample consisted of homogeneously-sized small grains, probably 
because of the increased surface roughness. Conversely, if a srniple of intermediate-size particles 
containeil a high lret|uency of tiny grains, then the overall threshold speed may be anomalously 
high. Thus care was taken to ensure relatively accurate size distributions for the test sample; 
nonetheless, it was not always possible to obtain well-graded samples. 
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Figure 12.- Comparison of the shape of walnut shells with losses and water-deposited silt. The shape, degree of 
sorting, and cleanliness of the walnut shells is very similar lo windblown losses. 


Test bed length and threshold determination ■ An important parameter in thresholil tests is the 
length of the bed of material being tested (fig. 13). If the bed length was not sufficiently long, the 
experimental value of the threshold speed would be unreulistically high. If the bed length was longer 
than a critical length, the experimentally determined value of 1’*^ became a constant, rurthermore. 
the critical minimum bed length is tiot constant for all particle diameters but appears to be a 
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Figure 13. Kffect of length of particle lest bed on threshold at I bar pressure for 74-l4‘>jum silica sand. (x;ngths 
shorter than about I .S in result in anomalously high friction thresliohl. 
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Iiinctioii ol ihc particle diameter and atmospljeric pressure, I'or tlie smaller-sized particles the 
reijiiired critical hed length is suhstantially shorter than that for the larger-sized particles. I'or the 
major size ranges of test particles, a series of tests was run to determine tlie critical minimum hed 
lengMi for the subsetiuent threshold tests. Typically a critical bed length of 2 m is needed for 
particles with diameters ot .SOOjum and about l-m bed length for smaller particles. 

Saltation threshold is defined here as occurring whui grains of a given size are set into motion 
over an entire bed. High-speed motion pictures and direct observations showed that typically in the 
wind tunnel, as the wind is slowly increased, a speed is reached at which a few grains begin to 
quiver, roll, and saltatc; this is followed by a Hurry of activity, then by relatively no motion. This is 
not considered true threshold since it appears to be the result of the manner of emplacement of the 
particles on the test bed. Evidently there are some particles that are “perched” above the general 
surface and are highly susceptible to movement by the wind. Thus the procedure in the tunnel is to 
raise the wind speed to the point that these grains are removed, with the surface of the bed being 
repositioned naturally by the wind. The velocity is decreased below this initial movement then 
slowly increased until true saltation ot the bed is observed. Several data points were taken at each 
pressure and several particle samples were tested two or more times in subsequent runs to determine 
repeatability of the data; the error band was found to be on the order of 20 percent. 

Saltation can be detected by several methods. Previous threshold experiments (refs. 14. 18, 
and 1 9) conducted under one-atmosphere pressure involved visual observations of the entire test 
bed (saltation results in a cloud of particles moving over the entire surface) and observation of a 
■^2-cm-sized patch ol the bed through a telescope. This procedure, however, is not practical in 
MARSWIT for low pressure experiments since the tunnel test bed is too far removed from the 
observation port of the control room (fig. 7). This saltation is detected by means of. first, a high 
resolution closed-circuit television system having a lOOO-mm focal length camera lens to observe a 
patch ot the test bed about 3 cm across; and second, a luser beam-photometer system (fig. 14) in 
which a laser beam is directed longitudinally down the test bed. reflected off a mirror at the end of 
the tunnel and directed back across the test bed where it activates a photocell. As soon as saltation 
occurs, the particles in motion interfere with the laser beam, causing less light to reach the photocell 


TUDE TO PROTECT 
PHOTOCELi FROM 
EXTRANEOUS LIGHT 


RETURN LASER 
BEAM MODULATED 
BY SALTATING 
PARTICLES 


AERODYNAMIC 

SHIELD 



l igiire 14. .Sek-miitic diagram showing iaser-bcam sallalion threshold device. A laser beam is directed across the 
floor ot the tunnel, immediately above the test particles, and rellected frtmi a mirror onto a photocell. As soon as 
particles begin saltation, they intertere with the laser beam which is detected by the photocell. 


aiui reiliicinp tlie signal generated by the photocell. The third method involves the >ise of an 
electrometer. It has been found in the wind tunnel that beds of saltating particles build up 
substantial electrical charges. An electrometer probe placed at the end of the test beil shows a rapid 
buililup of charge as soon as saltating particles begin tr> impinge upon the probe. Tiiis multiple 
method of saltation detection (television observation, laser beam-photometer, and electrometer) is 
used for both atmospheric and tow piessure tests for uniformity in iletection. 

Threshold as a I'unction of Reynolds Number 

Several recent calculations of saltation parameters (primarily threshold friction speeds) for 
Mars have been based on Bagnold's (ref. 21 ) terrestrial work. Many of these estimates (rels. 27 -29, 
and others) are based on the assumption that Bagnold’s coefficient /I is a unique function of the 
particle friction Reynolds number li (i.e.. A = AUi)X U A were a “universar’ function of B, then 
one curve should suffice for all conditions, including various atmospheric pressures. This would 
allow fairly direct extrapolations from experiments conducted at one atmosphere to tow atmo- 
spheric pressures such as Mars. However, parameters such as interparticle forces resulting from 
cohesion by adsorbed water or electrostatic charges which may be functions of atmospheric 
pressure would severely complicate the extrapolation. 

in a recent report, Sagan and Bagnold (ret. 30) argue that since Mars is comparatively dry, 
interparticle forces resulting from cohesion by absorbed water would be negligible. They assume 
that A is a certain function of Reynolds number B, ’out their I unction is not derived trom Bagnold s 
original experimental curve. Instead, they present a curve for air extrapolated from the experimen- 
tal threshold data of White (ref. 31) for particles in water. Since it is presumed that particles 
immersed in water are cohesionless (at least, as might be related to adsorbed moisture), their 
threshold curve for air is an estimate for cohesionless particles and deviates considerably from 
Bagnold's original curve for particles in air of diameter less than 100 /on. They present a threshold 
friction velocity curve for Mars in which there is no upturn in friction velocities for the smaller 
particles (ref. 30). However, with the return of Viking Lander pictures (fig. 3) showing rather 
cohesive, fine grained materials. Sagan agrees (personal communication) that some lorm ot inter- 
particle force must exist for small particles oi> Mars. In addition, as pointed out by Pollack vt al. 
(ref. 17). the fricl ionless case of Sagan and Bagnold would predict particle motion to be much more 
frequent than is known to be the case from Viking results. 

To determine the effect of interparticle force on threshold. Iversen ct al. (ref. conducted 
one-atmosphere wind-tunnel experiments for particles with specitic gravities ranging Irom 0.21 to 
11.35 and diameters ranging from 8 gm to I29()ium. giving a wide range of particle friction 
Reynolds numbers. These results were combined with the limited threshold data obtained at low 
pressure by Weinberger and Adlon (ref. 32). and it was found that the "universality" o1 A(B) lor 
small particles at low pressure was essentially invalid. 


i:\periinental Results Obtained Imder l ow Atmospheric Pressure 

Initial results for friction threshold speeds conducted under low atmospheric pressure are 
presented in reference 15. Ihe results presented here are expansions of those results and refine- 
ments to take into account such factors as the influence of bed length on threshold. 
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TABLE 1, MATERIALS USED IN LOW 
PRESSURE THRESHOLD EXPERIMENTS 

Material 

Size, pm 

Specific gravity 

Caldtc 

6 200 

2.7 

Copper 

12 60 

8.9 

Cork 

175 850 

0.27 

Olass microspheres 

5 710 

2.47 

Instant tea 

400 - 1200 

0.21 

Polystyrene DVB 

175 850 

1.01 

Sand 

20 - 1400 

2.67 

Silt 

<5 60 

2.65 

Walnut shell 

<20 1400 

1,34 




For eacli sample material tested 


performed under pressures ranging from 
4 mb (400 Pa) to 1000 mb (10* Pa), Fig- 
ure 1 5 shows the threshold velocity as a 


function of atmospheric pressure for 


three samples of walnut shells, sizes 
64 pm, 330 pm, and 630 pm. Similar 
graphs were obtained tor each sample 
material. Values of were taken from 
these graphs for each particle diameter 
to derive vs Dp curves for several 
given atmospheric pressures (fig, 16). 


1000 r 


100 h 


A curve for demonstrating the func- 
tional relationship between A and B and 
its “nonuniversaiity” at different pres- 
sures is the dimensionless friction speed: 
V*(liplpp- P)8>'y'^ vs 

dimensionless diameter: 

Dp[(pp-p)8lpv^]^'^ 

(ref. 19). Figure 17 is such a plot for the 
ground walnut shells. The existence of 
separate curves for each of the diameters 
indicates the presence of additional 
terms (i.e., interparticle forces, etc.) in 
the functional relations and confirms the 
existence of the “upturn” in threshold 
curve for Mars with the corresponding 
existence of an “optimum” grain size for 
Mars. 


10 h 


[i #21 WALNUT AVESIZE 630 urn 
#28 WALNUT AVESIZE 330 
A #35 WALNUT AVE SIZE 64 #m 
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FRICTION THRESHOLD SPEED, V.,. Cm/6C 


Piguie 1 5. - Curve showing increase in friction threshold 
speed (F« ) with decrease in pressure. 


Not all of the parameters (e.g., vis- 
cosities and temperatures) involved in 
aeoKan movement can be satisfied in 
wind-tunnel simulations conducted on 
Earth, even under low-pressure condi- 
tions. Thus a combination of theory and 
wind-tunnel results must be employed 
for extrapolation to Mars as discussed in 
reference 15. Figure i8 shows threshold 
results performed at low pressure but 
with Earth air; figure 16 is an extrapola- 
tion of the wind-tunnel results to Mars 
based on the expressions derived by 
Iversen ct al. (ref, 19). using the 
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Figure 16.- Martian particle threshold curves as a function of particle size at three pressures. Scale on left is in 
m/sec, two scales on right are equivalent free-stream (above the boundary layer) wind velocities ( Vgf), based on 
atmospheric models by Pollack et al. (refs. 16 and 17). Case I is for winds blowing over a flat smooth surface 
composed of erodible grains; Case 2 is for a surface containing cobbles and small boulders. 





Figure 17.- Dimensionless threshold friction speed as a function of dimensionless particle diameter for several 
different ground walnut shell distribution. Experiments were performed from a minimum pressure of S mb to 
1 atm. If Bagnold's A is uniquely only a function of 1) for nonterrestrial conditions, then only a single curve 
should suffice to fit all the data. Obviously, this is nut the case, as each separate particle diameter range exhibits 
its own curve and A = A{it) is not true when the pressure is varied. However, at a constant pressure, (/I ^/i>* is a 
unique function of (/i//l)’ . Tlie differences that occur owing to pressure changes may possibly be explained by 
clianging of electrostatic and interparticlc cohesive forces. 



Figure 18- Threshold frictlorr velocities of 
particle diameter for walnut shells (circles) in 
Earth air at low pressure, compared to 
theoretical threshold curves, calculated using 
Earth air and Earth gravity to simulate the 
conditions. 


appropriate values for kiiietnatie viscosity atul tlcn- 
sity for the martian atmosphere. 

Minimum threshold friction speed to 

move particles on Mars is about 2..S m/sec. Depend- 
ing upon the atmospheric conditions (ref. 17), at a 
nominal tnartian surface pressure of 5 mb (500 Pa) 
this would correspond to a wind velocity a couple 
of meters above the surface of about 125 m/sec for 
a smooth tlat surface without boulders, or to about 
43 m/sec for a surface with boulders (ref. 1 5). 

The particle size most easily moved on Mars 
by the wind is about 1 60 jum in diameter as shown 
in figure 1 6. In an active wind regime over an area 
of relatively fixed particle formation, this size 
particle (fine sand) would be easily mobilized and 
eventually removed, either leaving behind both 
finer (silt) and coarser materials, or possibly caus- 
ing the silt to go into suspension, depending on the 
ratio of sand to silt and their relative placement 
with respect to the wind. 


OTHER PARAMETERS AFFECTING AEOLIAN PROCESSES ON MARS 


The wind speeds needed to raise particles on Mars discussed under “Threshold Friction Speed 
Experiments” apply to relatively steady winds blowing across tlat smooth plains consisting of loose, 
uniformly size-graded particles. The series of experiments to obtain for these parametrically 
simple but probably unrealistic conditions, was a first-order analysis in which the number of factors 
was kept to a minimum. It is obvious from the Viking Lander pictures that for at least two localities 
on Mars, the surface is very complex; Viking Orbiter pictures indicate that the rest of the planet is 
equally or more complex. 


In this section we will discuss some of the complicating factors that would alter the values of 
wind speeds needed to raise particles beyond those described above. In general, these effects lower 
the tree-stream wind speeds needed to set particles into motion; however, very little research has 
been conducted on these effects, particularly in regard to low atmospheric pressure conditions, and 
the results presented here are exploratory in nature and should be taken as preliminary. 


Swiilhm Motion TiiivsUoUl 

Many authors (c.g., rcl . 33 ami others) have rccogni/.eil the rliflievilty of pmih.diig stirikieiUly 
high tree-stream wintts to raise particles on Mars. They have suggested that * dust-devil oi smular 
cyclonic-type winds may i>hiy an important or even domhtant role in imttatmg martian dust storms. 

In order to assess the relative elTectiveness of swirling moti(»n winds as a means of 
particle threshold, a series of exploratory experiments was cimducted at one 
determine the strength and size of vortices needed to hit erodihle particles (table >) 
surface. Research involving a vortex generator (ref. 34) has 

University Aerospace lingineering Department tor the past several years, The vo ex, tormeit w m 
its axis perpendicular to the ground or surface plane, is produced by passing air through a rotating 
le^^crb and stationary nLzle located some distance above the surface. A vortex can be formed 
by passage of air either upwards or downwards through the honeycomb or with no air passage at a 1. 
The vorttx -nail-scale model of atmospheric vortices such as the tornado 

or dust-devil (fig. 19), 

The swirling motion due to the atmospheric vortex (dustKlevil) causes a radial pressure 
gradient at the surface 


^ = - — - ~ 
dr r o 


^ Ar ^ Ap = pro^<^o^ 




where p is pressure, r is radial distance 
from the vortex center, p is air density, 
uq is tangential (swirl) speed, to,, is the 
angular speed of the generator honey- 
comb, and the characteristic radius r,, is 
as defined in equations (4) where A/) is 
the maximum surface pressure dillerence 
from the center of the vortex to that at 
large radius. 

The maximum pressure difference, 
A(), was measured on the surface plane 
of the vortex generator with a differen- 
tial pressure transducer for the range of 
angular speeds of the generator. The 


TABLt 2. MATERIALS USED IN WIND-TUNNEL 
TESTS _____ 


Material 

Density. ~gin/cnr | 

Diameter. ~Min 

O Instant tea 

0.21 

719 

^ Silica gel 

0.89 

17; 189 

□ Nut shell 

1.1 

40 to 359 

^ Clitvor seed 

1.3 

1290 

"V* Sugar 

1.59 

393 

^ Glass 

2.42 

31 to 48 

O Glass 

, 2.5 

38 ti) 586 

L Sand 

2.65 

526 

Q Aluiniimiti 

2.7 

36 to 204 

A Glass 

3.9') 

55 to 519 

O Cupper oxide 

6.0 

10 

V Hrt)ii/e 

7.8 

1 

6l6 

() Cupper 

8.‘»4 

12; 37 

O l.ead 

11.35 
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ORIGINAJ, I'AfilO KS 
OP POOH oi Ai.m- 

values o!' clKuacterisiic radius calcu- 
lated frojji t!ie ineasureil values of A/i are 
shown in nj?ure 2(1. A linear lit ol the 
data results in an c(|uation lor 

-- 0.0002 ^ 6.H7 (5) 

where is in cm ami in rad/sec. 

Assume tliat the lop layer of par- 
tides of thickness ei|ual to particle diam- 
eter Dp is lifted by the vortex at thresh- 
old. 1 he forces actint* on this layer per 
unit urea are then us shown 


In the figure above, Ap represents the 
instantaneous pressure difference 
between the lower and upper surfaces of 
the particle layer, caused by the vortex 

Flsure 19... Vtew »f Ihe »i,line (cycUmlc) mollon d.vi«, P“' '>'5 where 

raising a small “duM dovir of parllclaa above a leal bad. ore obsom-d to lift ofl tire nut- 

face. The factor a is included to account 

for pressure relief that might occur prior to particle lift-off; oj is the stress due to interparticle 
force, T is an upward normal stress due to the effects of viscosity, ppgDp is the weight of the 
particle layer per unit area, and ct is an unknown constant of proportionality. If the upper layer 
were completely sealed a would be 1; if there were leakage « would be less than 1. It the surtace 
stress T is defined as t = , then the value of is threshold where forces are in 

equilibrium can be written as 


'■o^o 


(«4 






The characteristic threshold speed is not a true speed but is proportional to the square root of 
the pressure difference Ap (i.v.. ~ ( A/>/p)' )• 

Unfortunately neither the function of^Dp) nor the function nor the ratio a are 

known, and sufficient data to delennine these functions empirically is not yet available. The values 
of at threshold are shown as a function of particle diameter in figure 21 and as a function of 
(p.xDJp)' in figure 22. |■igure 22 does show a trend for increasing angular threshold speed for 



Figure 20. C'luiractcfistic radius of Ore vortex as g I 
a I'uttctioii of •’otor speed (see cqs. (4) and (S)). o 4 
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Figure 21.- Characteristic vortex thresh- 
hold speed as a function of par- 
ticle diameter. 
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Figure 22.- Characteristic vortex thresh- 
old speed as a func- 

tion of particle parameter. Note that 
for large particles the pressure differ- 
ence Ap Is neatly equal to the weight 
per unit area of a single layer of par- 
ticles. For small particles the pressure 
difference ranges up to nearly ten 
times the weight per unit area of a 
single layer. One explanation for the 
difference is small particle cohesion 
(intcrparticlc stress o»). 
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Iiirtid viilufs of this function. However, since the effects of Reynoltls miinbei and 

intei'particle force are not accounted for. there is a lar^e amount of data scatter for most of the 
range of the data. 

The entrainmem mechanism fm raising surface particles with a vertical vortex (dust devil) is 
obviously quite different from the uniform wind boundary-layer case. The variation with particle 
diameter is much less pronouneed for vorte.x threshold than it is for boundary-layer IhreshoUl, I'he 
important points to note from these exploratory experiments is that the si/e and density seem to 
have little importance in swirling motion threshold except for the two extreme eases of large lead 
ami bron/e particles. Thus cyclonic motion winds would appear to be very effective in raising 
erodible particles of all types and si/.cs. Once lifted above the surface, much lower winds than 
“normal" threshold could continue to carry particles along. 


Behavior of Windblown Particles at I.ow Pressure 


Once particles are set into saltation by the wind, their behavior at low pressure is quite 
different from that of particles under one atmosphere conditions. !'or example, because much 
stronger winds are required for threshold, the particles carried by the wind in saltation and 
suspension will also be moving faster after threshold. When they impact other particles and rock 
surfaces they will impart more energy to the impacted surface, which may result in a greater erosive 
capability at low pressure as well as potentially lowering the dynamic threshold of the impacted 
grains. 

Moreover, partly because of the increased particle velocities and reduced atmospheric drag, the 
saltation path length has been calculated tt> be up to 50 times longer than on lianh. according to 
the assumed surface conditions. These values arc based on numerical simulations, using wind-tunnel 
data obtained at one atmosphere surface pressure (ref. 5). 


Since the geometry of ripple length, erosional scour patterns, and other aeolian features is 
dependent on the characteristics of the saltation tiajectories. we should expect to see differences in 
these features on Mars. I'or example. Viking Lander photographs show a pronounced scour /.one in 
the deposits of fine-grained particles around many of the rocks: although similar patterns develop 
around rocks on liarth. there appear to be geometric differences in the martian cases, primarily in 
the I'orm of deeiu’r and longer scour /ones. This may be explained by the differences in the 
boundary layers near the surface. The properties of the viscous sublayer flow within the boundary 
layer influences the motion of the grains. On Larth. the typical range of saltation occurs when there 
is no viscous sublayer or at most a very minute one t'-l mm thick; ref. 5). The viscous subhr/er is 
proportional to which can be iqi to lO times as large as that on liarth or ~l cm thick. I hus the 
effect is to have greater viscous interactions occurring at the surface on Mars and in the low' pressure 
experiments, and the results are deeper erosion and extension of erosion further downstream, 


l o test potential causes of such differences, several experiments were conducted in MARSWI T. 
A rock about I 2 cm across was placed in the tunnel on a bed of fine sand; under one atmosphere 
condition, the sand was eroded to determine (he sct>ur /ones in relation to the rock. Vortices and 
turbulence around the rock resulted in eroded /ones in front of the rock and oti its leeward sides, a 
pattern tt> be expected from the horseshoe vortex flow lield described (refs. 4 and 1 4). I he same 
test was run under a surface pressure of 5..^ mb. 
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AlthmiMli thf h;isic p:it!om Wiis tiu! s;iim- lor both LiisfN(hoih result liitin liorsesline voitievs) 
the erosion is more proiunmeeil lor Ihe siimilateil martian ease. Holh eases *vere run al wimi speeds 
slittlilly hitthei Hum threshold; heeaiise is many times hittlier lor nuiilian eonililioiis than lor 
I'.arth eoiiditions. Hie rree-stivam wind veloeity is also eonsiderahly hir'Jier. 1 he pressure diHeienlial 
iietintt on the partieles is approsimately the same, however. Onee Ihe p.utieles are jiieked up hy Hie 
wind, they also travel I'asler. Wlien the saltatinn grains strike ilie loek. they bounee oil the laeeol 
the rock with pivaler energies and erode the /one in Ironl of the roek to a greater denree. Moreover, 
because the saltation path lenuth is {tivater under martian conditions, as the plains are swept into 
the horseshoe vortex around the rock, tliey will strike the snrlaee at a prealer distance downwind 
Irom the rock than they did at one atmosphere. Such a result has the elfect of leiiptheninp tin 
aeohan "shadow” zone behind rocks and other obstructions. Since the plains are travelinp al a 
prealer velocity and liavc a larper momentum, upstream plains are also less likely to be "captured" 
hy the rock as they enter its flow t'ield. 

Anotlier aspect of particle behavior at low pressure that affeets thresliold and may alTecl rates 
ol erosion also, is the spin rate of grains in saltation. Although it has been known for some time that 
saltatmg grains are spinning, the significance of that fact has not been fully appreciated, l or 
example. White and Schulz (ref. 20) have shown that part of the lift on the particle in the saltation 
trajectory is attributable to the Magnus effect resulting from spinning grains. As the particles begin 
to htt above the bed, they begin spinning due to shear forces in the boundary layer. Analysis of 
high-speed motion pictures of saltating grains show that they spin up to 200 rps under one 
atmosphere condition. Hxtrapolation to martian conditions show that spin rates are more than 
doubled, or about 450 rps. Although the effect of rapidly spinning grains on rates of erosion has not 
oeen assessed, it may be significant. Future experiments will utilize MARSWIT to test the rate of 
spin at low pressure as well as to assess the effect of spinning grains on rales of erosion. 

SUMMARY AND CONCLUSIONS 


That aeolian processes have played and continue to play a significant role in the evolution of 
the surtace of Mars cannot be disputed. Unfortunately, aeolian processes on Mars are pooriy 
understood, primarily because ot the drastically dilferent aeolian environment compared to I'arth's; 
yet, without good knowledge ol these processes, the geological history of Mars remains shadowed 
Rates of aeolian erosion provide a case in point. Prior to the Viking mission, orbital pictures from' 
Mariner 9 showed numerous features that were attributed to winds, including dunes, crater streaks, 
and other depositional features. Pedestal craters, yardangs. sculpted laminated terrain and other 
features observed from orbit prompted the general conclusion that vast regions had been eroded to 
significant depths by the wind. It was reasoned that the frequent dust storms and the postulated 
high winds would be capable of accomplishing high rates of erosion. 

Viking lesults, however, are raising serious questions about just how much wind erosion lias 
occurred. Pictures from orbit, for example, show that many of the pedestal craters are not erosional 
but are nrimary lamlforms some are fresh impact craters tref. .L5) and others are probably 
volcanoes, although som still at'peai' to result from erosion. 

Viking Lander pictures in ('hryse Planitia lead to the conclusion that the plains may have 
undergone from 1 10 in of surlace erosion since their emplacement, a rather small amount 



coiisiilcrinp tlicir yiviit at?o as tentatively determined from crater counts (Binder, personal communi- 
cation). Moreover, rocks visible in tlic Lamler pictures show few of the surface textures and 
features expected from wind erosion, as least hy terrestrial comparisons. Similar conclusions have 
been reached for the Vikinp II Lander site. 

Other surprises from Vikinti are causiiiR the questions about wind erosion on Mars to be 
reasked The basic question is: How effective is wind erosion on Mars under the present environ- 
ment and how might it have been different in the past? To answer this question, one must know: 


1 . The range of wind speeds at which particles are driven by the wind. 

2 . The characteristics of saltation in the martian atmosphere. 

3. The effects of electrostatically-charged particles. 

4. Rock weathering in the martian environment. 

5. Frequency-magnitude determinations ot martian surface winds. 


Solutions to these eoiuplex problems will undoubtedly take yearn to derive. Nonetheless 
results presented here are focusing on these problems, which, when combined with analysis of 
spacecraft data will lead to a much clearer understanding of aeolian processes on Mars. 
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